Neutron Star Interiors and the Equation of 
State of Superdense Matter 



Fridolin Weber 1 , Rodrigo Negreiros 2 , and Philip Rosenfield 3 

1 San Diego State University fweber@sciences.sdsu.edu 

San Diego State University negreiro@sciences.sdsu.edu 
3 San Diego State University philrose@sciences.sdsu.edu 

Summary. Neutron stars contain matter in one of the densest forms found in the 
Universe. This feature, together with the unprecedented progress in observational 
astrophysics, makes such stars superb astrophysical laboratories for a broad range 
of exciting physical studies. This paper gives an overview of the phases of dense 
matter predicted to make their appearance in the cores of neutron stars. Particular 
emphasis is put on the role of strangeness. Net strangeness is carried by hyperons, 
K-mesons, H-dibaryons, and strange quark matter, and may leave its mark in the 
masses, radii, moment of inertia, dragging of local inertial frames, cooling behavior, 
surface composition, and the spin evolution of neutron stars. These observables play 
a key role for the exploration of the phase diagram of dense nuclear matter at high 
baryon number density but low temperature, which is not accessible to relativistic 
heavy ion collision experiments. 



1 Introduction 

Neutron stars are dense, neutron-packed remnants of stars that blew apart in 
supernova explosions. Many neutron stars form radio pulsars, emitting radio 
waves that appear from the Earth to pulse on and off like a lighthouse beacon 
as the star rotates at very high speeds. Neutron stars in x-ray binaries ac- 
crete material from a companion star and flare to life with a burst of x-rays. 
The most rapidly rotating, currently known neutron star is pulsar PSR J1748- 
2446ad, which rotates at a period of 1.39 ms (which corresponds to a rotational 
frequency of 719 Hz) [1]. It is followed by PSRs B1937+21 [2] and B1957+20 
[3] whose rotational periods are 1.58 ms (633 Hz) and 1.61 ms (621 Hz), re- 
spectively. Finally, the recent discovery of X-ray burst oscillations from the 
neutron star X-ray transient XTE J1739-285 [4] could suggest that XTE 
J 1739-285 contains the most rapidly rotating neutron star yet discovered. 
Measurements of radio pulsars and neutron stars in x-ray binaries comprise 
most of the neutron star observations. Improved data on isolated neutron 
stars (e.g. RX J1856. 5-3754, PSR 0205+6449) are now becoming available, 
and future investigations at gravitational wave observatories focus on neutron 
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Fig. 1. Neutron star compositions predicted by theory. 

stars as major potential sources of gravitational waves (see [5] for a recent 
overview). Depending on star mass and rotational frequency, the matter in 
the core regions of neutron stars may be compressed to densities that are up to 
an order of magnitude greater than the density of ordinary atomic nuclei. This 
extreme compression provides a high-pressure environment in which numerous 
subatomic particle processes are likely to compete with each other [6, 7]. The 
most spectacular ones stretch from the generation of hyperons and baryon 
resonances (S, A, S, A), to quark (u, d, s) deconfmement, to the formation of 
boson condensates (ir~ , K~, H-matter) [6, 7, 8, 9, 11, 12] (see Fig. 1). In 
the framework of the strange matter hypothesis [15, 16, 17], it has also been 
suggested that 3-flavor strange quark matter-made of absolutely stable u, d, 
and s quarks-may be more stable than ordinary atomic nuclei. In the latter 
event, neutron stars should in fact be made of such matter rather than ordi- 
nary (confined) hadronic matter [18, 19, 20]. Another striking implication of 
the strange matter hypothesis is the possible existence of a new class of white- 
dwarfs-like strange stars (strange dwarfs) [21]. The quark matter in neutron 
stars, strange stars, or strange dwarfs ought to be in a color superconduct- 
ing state [22, 23, 24, 25]. This fascinating possibility has renewed tremendous 
interest in the physics of neutron stars and the physics and astrophysics of 
(strange) quark matter [11, 22, 23]. This paper discusses the possible phases 
of ultra-dense nuclear matter expected to exist deep inside neutron stars and 
reviews selected models derived for the equation of state (eos) of such matter 
(see also [6, 7, 8, 9, 10, 11, 12, 13, 14]). 




Fig. 2. Models for the equation of state (pressure versus energy density) of neu- 
tron star matter [11]. The notation is as follows: RMF=relativistic mean-field 
model; DD-RBHF=density dependent relativistic Brueckner-Hartree-Fock model; 
n=neutrons; p=protons; H=hyperons, K=K~ [u, s] meson condensate; Q=u, d, s 
quarks; H-matter=H-dibaryon condensate. 

2 Neutron Star Masses 

In 1939, Tolman, Oppenheimer and Volkoff performed the first neutron star 
calculations, assuming that such objects are entirely made of a gas of non- 
interacting relativistic neutrons [26, 27]. The eos of such a gas is extremely 
soft, i.e. very little additional pressure is gained with increasing density, as 
can be seen from Fig. 2, and predicts a maximum neutron star mass of just 
0.7 Mq (Fig. 3) at an unrealistically high density of 17 times the density 
of nuclear matter (Fig. 4). It is interesting to note that the inclusion of 
interactions among the neutrons increases the star's maximum mass from 
0.7 Mq to around 3 Mq (Figs. 3 and 4). However, the radii of the latter 
are so big that mass shedding from the star's equator occurs at rotational 
frequencies that are considerably smaller than those observed for PSR J1748- 
2446ad, 716 Hz (1.39 ms) [1], or B1937+21, 630 Hz (1.58 ms) [2]. An in- 
teracting neutron gas thus fails to accommodate the observed rapidly ro- 
tating neutron stars. The other extreme, a non-interacting relativistic neu- 
tron gas, fails too since it does not accommodate the Hulse- Taylor pul- 
sar (M — 1.44 Mq) [28], and also conflicts with the average neutron star 
mass of 1.350 ± 0.004 Mq derived by Thorsett and Chakrabarty [29] from 
observations of radio pulsar systems. More than that, recent observations 
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indicate that neutron star masses may be as high as around 2 Mq. Ex- 
amples of such very heavy neutron stars are Mj 75i+i807 = 2.1 ± 0.2 Mq 
[30], M 4U16 36 +5 36 = 2.0 ±0.1 Mq [31], Mveiax-i = 1.86 ± 0.16 M Q [32], 
Mcygx-2 = 1-78 ± 0.23 Mq [33, 34]. Large masses have also been reported 
for the high-mass x-ray binary 4U 1700-37 and the compact object in the 
low-mass x-ray binary 2S0921-630, M 4U 1700-37 = 2.44 ± 0.27 M Q [35] and 
M2S0921-630 = 2.0 — 4.3M Q [36]. respectively. The latter two objects may 
be either massive neutron stars or low-mass black holes with masses slightly 
higher than the maximum possible neutron star mass of <~ 3M Q . This value 
follows from a general, theoretical estimate of the maximal possible mass of 
a stable neutron star as performed by Rhoades and Ruffini [37] on the basis 
that (1) Einstein's theory of general relativity is the correct theory of grav- 
ity, (2) the eos satisfies both the microscopic stability condition dP/de > 
and the causality condition dP/de < c 2 , and (3) that the eos below some 
matching density is known. From these assumptions, it follows that the maxi- 
mum mass of the equilibrium configuration of a neutron star cannot be larger 
than 3.2 Mq. This value increases to about 5M if one abandons the causal- 
ity constraint dP/de < c 2 [38, 39], since it allows the eos to behave stiffer at 
asymptotically high nuclear densities. If cither one of the two objects 4U 1700- 
37 or 2S0921-630 were a black hole, it would confirm the prediction of the 
existence of low- mass black holes [40]. Conversely, if these objects were mas- 
sive neutron stars, their high masses would severely constrain the eos of dense 
nuclear matter. 



3 Composition of Cold and Dense Neutron Star Matter 

A vast number of models for the equation of state of neutron star matter has 
been derived in the literature over the years. These models can roughly be 
classified as follows: 

• Thomas- Fermi based models [41, 42] 

• Schroedinger-based models (e.g. variational approach, Monte Carlo tech- 
niques, hole line expansion (Brucckner theory), coupled cluster method, 
Green function method) [8, 43, 44, 45, 46, 47, 48] 

• Relativistic field-theoretical treatments (relativistic mean field (RMF), 
Hartree-Fock (RHF), standard Brueckner-Hartree-Fock (RBHF), density 
dependent RBHF (DD-RBHF) [49, 50, 51, 52, 53, 54] 

• Nambu-Jona-Lasinio (NJL) models [55, 56, 57, 58, 59, 60] 

• Chiral SU(3) quark mean field model [61]. 

A collection of equations of state computed for several of these models is shown 
in Fig. 2. Mass-radius relationships of neutron stars based on these equations 
of state are shown in Fig. 3. Any acceptable nuclear many-body calculation 
must correctly reproduce the bulk properties of nuclear matter at saturation 
density, no = 0.16 fm~ 3 . These are the binding energy, E/A = —16.0 McV, 
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Fig. 3. Mass- radius relationship of neutron stars and strange stars [11]. The strange 
stars may be enveloped in a crust of ordinary nuclear material whose density is below 
neutron drip density [18, 62, 63]. 

effective nucleon mass, = 0.79 ton, incompressibility, K ~ 240 MeV, and 
the symmetry energy, a s — 32.5 MeV. 

3.1 Hyperons and baryon resonances 

At the densities in the interior of neutron stars, the neutron chemical potential, 
/i", is likely to exceed the masses, modified by interactions, of £, A and possi- 
bly S hyperons [64]. Hence, in addition to nucleons, neutron star matter may 
be expected to contain significant populations of strangeness carrying hyper- 
ons. The thresholds of the lightest baryon resonances (A~ , A , A + , A ++ ) are 
not reached in relativistic mean- field (Hartree) calculations. This is different 
for many-body calculations performed at the relativistic Brueckner-Hartree- 
Fock level, where A's appear very abundantly [65]. In any event, pure neutron 
matter constitutes an excited state relative to hyperonic matter which, there- 
fore, would quickly transform via weak reactions like 



to the lower energy state. The chemical potentials associated with reaction 
(1) in equilibrium obey the relation 



where py" = since the mean free path of (anti) neutrinos is much smaller 
than the radius of neutron stars. Hence (anti) neutrinos do not accumulate 
inside neutron stars. This is different for hot proto-neutron stars [66]. Equa- 
tion (2) is a special case of the general relation 



n —>■ p + e +v e 



(1) 




(2) 
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Fig. 4. Neutron star mass versus central density (in units of nuclear matter density, 
e = 140 MeV/fm 3 ). 



fi x = B x n n - q x l i e , (3) 

which holds in any system characterized by two conserved charges. These are 
in the case of neutron star matter electric charge, q x , and baryon number 
charge, B x . Application of Eq. (3) to the A hypcron (B A = 1, q A = 0), for 
instance, leads to 

/j, A = f/ 1 . (4) 

Ignoring particle interactions, the chemical potential of a relativistic particle 
of type x is given by 

fi x = uj(k Fx ) = + k 2 Fx , (5) 

where uj(kp x ) is the single-particle energy of the particle and kp x its Fermi 
momentum. Substituting (5) into (4) leads to 

kF 3 

kF n > \/ m A ~ m n — 3 fm _1 => n = n ~ 6n , (6) 

» 37T 

where tua = 1116 MeV and m n — 939 MeV was used. That is, if interactions 
among the particles are ignored, neutrons are replaced with A's in neutron 
star matter at densities of around six times the density of nuclear matter. 
This value is reduced to about 2 no by the inclusion of particle interactions 
[64]. Densities that small are easily reached in the cores of neutron stars. 
Hence, in addition to nucleons and electrons, neutron stars may be expected to 
contain considerable populations of strangeness-carrying A hyperons, possibly 
accompanied by smaller populations of the charged states of the £ and E 
hyperons [64]. Depending on the star's mass, the total hyperon population 
can be very large [64], which is illustrated graphically in Figs. 5 and 6 for 




rotating neutron stars whose equation of state is computed in the framework 
of the relativistic DD-RBHF formalism [52] . Aside from chemical equilibrium, 
the condition of electric charge neutrality of neutron star matter, 

J2 <1 X 4 X + 3 ir 2 n M 0{p M -m M ) = 0, (7) 

X=P,S ± ,E~ ,A++ ,...;e _ 

where M stands for ir~ or K~ mesons, plays a key role for the particle com- 
position of neutron star matter too. The last term in (7) accounts for the 
possible existence of either a tt~ or a K~ meson condensate in neutron star 
matter, which will be discussed in more detail in Sect. 3.2. Before, however, we 
illustrate the importance of Eqs. (2) and (7) for the proton-neutron fraction 
of neutron star matter. The beta decay and electron capture processes among 
nucleons, n — > p + e~ + D e and p + e~ — > n + v e respectively, also known 
as nucleon direct Urea processes, are only possible in neutron star matter 
if the proton fraction exceeds a certain critical value [67]. Otherwise energy 
and momentum can not be conserved simultaneously for these reactions so 
that they are forbidden. For a neutron star made up of only nucleons and 
electrons, it is rather straightforward to show that the critical proton frac- 
tion is around 11%. This follows from kp n = kp p + hp a combined with the 
condition of electric charge neutrality of neutron star matter. The triangle 
inequality then requires for the magnitudes of the particle Fermi momenta 
kp n < kp p + kp e , and charge neutrality dictates that kp p = fc_F e . Substituting 
kp p = kp e into the triangle inequality leads to kp n < 2kp p so that for the 
particle number densities of neutrons and protons n n < &n p . Expressed as a 
fraction of the system's total baryon number density, n = n p + n n , one thus 
arrives at n p /n > 1/9 ~ 0.11, which is the figure quoted just above. Medium 
effects and interactions among the particles modify this value only slightly 
but the presence of muons raises it to about 0.15. Hyperons, which may exist 
in neutron star matter rather abundantly, produce neutrinos via direct Urea 
processes like E~ — * A + e~ + v e and A + e~ — * £~ + v e [68]. The direct 
Urea processes are of key importance for neutron star cooling (see D. Page's 
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contribution elsewhere in this volume). In most cases, the nucleon direct Urea 
process is more efficient than the ones involving hyperons [69, 70]. 

3.2 Meson condensation 

The condensation of negatively charged mesons in neutron star matter is 
favored because such mesons would replace electrons with very high Fermi 
momenta. Early estimates predicted the onset of a negatively charged pion 
condensate at around 2n (see, for instance, Ref. [71]). However, these esti- 
mates are very sensitive to the strength of the effective nucleon particle-hole 
repulsion in the isospin T = 1, spin S = 1 channel, described by the Landau 
Fermi-liquid parameter g' , which tends to suppress the condensation mecha- 
nism. Measurements in nuclei tend to indicate that the repulsion is too strong 
to permit condensation in nuclear matter [72, 73]. In the mid 1980s, it was 
discovered that the in-mcdium properties of K~ [us] mesons may be such that 
this meson rather than the ir~ meson may condense in neutron star matter 
[74, 75, 76]. 

The condensation is initiated by the schematic reaction e — ► K + v e . 
If this reaction becomes possible in neutron star matter, it is energetically 
advantageous to replace the fermionic electrons with the bosonic K~ mesons. 
Whether or not this happens depends on the behavior of the K~ mass, m* K - , 
in neutron star matter. Experiments which shed light on the properties of 
the K~ in nuclear matter have been performed with the Kaon Spectrometer 
(KaoS) and the FOPI detector at the heavy-ion synchrotron SIS at GSI [77, 
78, 79, 80, 81]. An analysis of the early K~ kinetic energy spectra extracted 
from Ni+Ni collisions [77] showed that the attraction from nuclear matter 
would bring the K~ mass down to m* K _ ~ 200 MeV at densities ~ 3n . For 
neutron-rich matter, the relation m* K _ /m K - ~ 1 — 0.2n/n was established 
[82, 83, 84], with mjf = 495 MeV the K~ vacuum mass. Values of around 
m* K _ ~ 200 MeV may be reached by the electron chemical potential, /x e , in 
neutron star matter [7, 64] so that the threshold condition for the onset of K~ 
condensation, \i e = m* K might be fulfilled for sufficiently dense neutron stars, 
provided other negatively charged particles {S~ , A~ , d and s quarks) arc not 
populated first and prevent the electron chemical potential from increasing 
with density. 

We also note that K~ condensation allows the conversion reaction n — ► 
p + K~ . By this conversion the nucleons in the cores of neutron stars can 
become half neutrons and half protons, which lowers the energy per baryon 
of the matter [85]. The relative isospin symmetric composition achieved in 
this way resembles the one of atomic nuclei, which are made up of roughly 
equal numbers of neutrons and protons. Neutron stars are therefore referred 
to, in this picture, as nucleon stars. The maximum mass of such stars has been 
calculated to be around 1.5 M [86]. Consequently, the collapsing core of a 
supernova, e.g. 1987A, if heavier than this value, should go into a black hole 
rather than forming a neutron star, as pointed out by Brown et al. [40, 82, 83]. 
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This would imply the existence of a large number of low-mass black holes 
in our galaxy [40]. Thielemann and Hashimoto [87] deduced from the total 
amount of ejected 56 Ni in supernova 1987A a neutron star mass range of 
1.43 — 1.52 Mq. If the maximum neutron star mass should indeed be in this 
mass range (~ 1.5 Mq), the existence of heavy neutron stars with masses 
around 2 Mq (Sect. 2) would be ruled out. Lastly we mention that meson 
condensates lead to neutrino luminosities which are considerably enhanced 
over those of normal neutron star matter. This would speed up neutron star 
cooling considerably [86, 70]. 

3.3 H-matter and exotic baryons 

A novel particle that could be of relevance for the composition of neutron star 
matter is the H-dibaryon (H=([itd] [ds] [su])), a doubly strange six-quark com- 
posite with spin and isospin zero, and baryon number two [88]. Since its first 
prediction in 1977, the H-dibaryon has been the subject of many theoretical 
and experimental studies as a possible candidate for a strongly bound exotic 
state. In neutron star matter, which may contain a significant fraction of A 
hyperons, the yl's could combine to form H-dibaryons, which could give way 
to the formation of H-dibaryon matter at densities somewhere above ~ 4 uq 
[89, 90, 91]. If formed in neutron stars, however, H-matter appears to unstable 
against compression which could trigger the conversion of neutron stars into 
hypothetical strange stars [90, 92, 93]. 

Another particle, referred to as exotic baryon, of potential relevance for 
neutron stars, could be the pentaquark, + ([ud] 2 s), with a predicted mass 
of 1540 MeV. The pentaquark, which carries baryon number one, is a hy- 
pothetical subatomic particle consisting of a group of four quarks and one 
anti-quark (compared to three quarks in normal baryons and two in mesons), 
bound by the strong color-spin correlation force (attraction between quarks 
in the color 3 C channel) that drives color superconductivity [94, 95] . The pen- 
taquark decays according to 6> + (1540) — ► K + [su] + n[udd] and thus has the 
same quantum numbers as the K + n. The associated reaction in chemically 
equilibrated matter would imply [i e = /j, K + /i n . 

3.4 Quark deconfinement 

It has been suggested already many decades ago [96, 97, 98, 99, 100, 101, 102, 
103] that the nucleons may melt under the enormous pressure that exists in the 
cores of neutron stars, creating a new state of matter know as quark matter. 
From simple geometrical considerations it follows that for a characteristic 
nucleon radius of ~ 1 fm, nucleons may begin to touch each other in nuclear 
matter at densities around (47rr^ r /3) _1 ~ 0.24 fm~ 3 = 1.5 no, which is less 
than twice the density of nuclear matter. This figure increases to <~ 11 no for 
a nucleon radius of rjv — 0.5 fm. One may thus speculate that the hadrons of 
neutron star matter begin to dissolve at densities somewhere between around 
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2 — 10 no, giving way to unconfined quarks. Depending on rotational frequency 
and neutron star mass, densities greater than two to three times n are easily 
reached in the cores of neutron stars so that the neutrons and protons in the 
cores of neutron stars may indeed be broken up into their quark constituents 
[6, 7, 11, 104]. More than that, since the mass of the strange quark is only m s ~ 
150 MeV, high-energetic up and down quarks will readily transform to strange 
quarks at about the same density at which up and down quark deconfinement 
sets in. Thus, if quark matter exists in the cores of neutron stars, it should 
be made of the three lightest quark flavors. A possible astrophysical signal 
of quark deconfinement in the cores of neutron stars was suggested in [105]. 
The remaining three quark flavors (charm, top, bottom) are way to massive 
to be created in neutron stars. For instance, the creation of charm quark 
requires a density greater than 10 17 g/cm 3 , which is around 100 times greater 
than the density reached in neutron stars. A stability analysis of stars with a 
charm quark population reveals that such objects are unstable against radial 
oscillations and, thus, can not exist stably in the universe [7, 11]. The same is 
true for ultra-compact stars with unconfined populations of top and bottom 
quarks, since the pulsation eigen-equations are of Sturm-Liouville type. 

The phase transition from confined hadronic (H) matter to deconfincd 
quark (Q) matter is characterized by the conservation of baryon charge and 
electric charge. The Gibbs condition for phase equilibrium then is that the 
two associated chemical potentials, pL n and /i e , and the pressure in the two 
phases be equal [6, 104], 

PH^ n ,^Ax},T) = P Q (^,^,T), (8) 

The quantity Ph denotes the pressure of hadronic matter computed for a 
given hadronic Lagrangian Cm({x}), where {%} denotes the field variables 
and Fermi momenta that characterize a solution to the field equations of 
confined hadronic matter, 

(i 7 ^ - m x )Vx(z) = J2 r Mx M(x)4> x (x), (9) 

M—(T,L0.7T.... 

(d^ + ml)a(x)= Yl r »xW*W. (10) 

X=p,n,S,... 

plus additional equations for the other meson fields (M = w, n, p, ...). The 
pressure of quark matter, Pq, is obtainable from the bag model [106, 107]. 
The quark chemical potentials /j, u , fi d , fi s are related to the baryon and charge 
chemical potentials as 

» u = \v n -\» e , V? = V* = \» n + \» e ■ (11) 

Equation (8) is to be supplemented with the two global relations for con- 
servation of baryon charge and electric charge within an unknown volume V 
containing A baryons. The first one is given by 



Neutron Star Interiors and the eos of Superdense Matter 



11 




3 6 9 12 15 18 3 6 9 12 15 18 3 6 9 12 15 
r (km) r (km) r (km) 

Fig. 7. Dependence of neutron star composition on spin frequency, v, for three 
sample compositions (left: hyperon composition, middle: quark-hybrid composition, 
right: quark- hybrid composition with quark matter in the color-flavor locked (CFL) 
phase [108]). The non-rotating stellar mass in each case is 1.4 Mq. vy. denotes the 
Kepler (mass-shedding) frequency, which sets an absolute limit on stable rotation. 



n=y = (l-r;)n H (^M e ,T)+ WM",M e ,T), (12) 

where rj = Vq/V denotes the volume proportion of quark matter, Vq, in 
the unknown volume V, and «h and uq are the baryon number densities of 
hadronic matter and quark matter. Global neutrality of electric charge within 
the volume V can be written as 

= | = (1 - V) MM", /A T) + V g Q (M", H e , T) + q L , (13) 

with qi the electric charge densities of hadrons, quarks, and leptons. For a 
given temperature, T, Eqs. (8) to (13) serve to determine the two independent 
chemical potentials and the volume V for a specified volume fraction r] of the 
quark phase in equilibrium with the hadronic phase. After completion Vq is 
obtained as Vq = r\V . Because of Eqs. (8) through (13) the chemical potentials 
depend on the proportion rj of the phases in equilibrium, and hence so also all 
properties that depend on them, i.e. the energy densities, baryon and charge 
densities of each phase, and the common pressure. For the mixed phase, the 
volume proportion of quark matter varies from < rj < 1 and the energy 
density is the linear combination of the two phases [6, 104], 

e = (1 - r?) e H ^", M e , M, T) + r, e Q ( M n , /A T) . (14) 

Hypothetical neutron star compositions computed along the lines described 
above are shown in Fig. 7. Possible astrophysical signals associated with quark 
deconfinemcnt, the most striking of which being "backbending" of isolated 
pulsars, are discussed in [6, 7, 11, 109, 110]. 
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3.5 Color-superconductivity 

There has been much recent progress in our understanding of quark matter, 
culminating in the discovery that if quark matter exists it ought to be in a 
color superconducting state [22, 23, 24, 25]. This is made possible by the strong 
interaction among the quarks which is very attractive in some channels. Pairs 
of quarks are thus expected to form Cooper pairs very readily. Since pairs of 
quarks cannot be color-neutral, the resulting condensate will break the local 
color symmetry and form what is called a color superconductor. The phase 
diagram of such matter is expected to be very complex [22, 23]. The complexity 
is caused by the fact that quarks come in three different colors, different 
flavors, and different masses. Moreover, bulk matter is neutral with respect to 
both electric and color charge, and is in chemical equilibrium under the weak 
interaction processes that turn one quark flavor into another. To illustrate 
the condensation pattern briefly, we note the following pairing ansatz for the 
quark condensate [111], 

mS^h) ~ A ^*f«hi + ^e af32 e fafb 2 + A 3 e^e fafb3 , (15) 

where V/ Q 1S a quark of color a = (r,g,b) and flavor f a = (u,d,s). The 
condensate is a Lorentz scalar, antisymmetric in Dirac indices, antisymmetric 
in color, and thus antisymmetric in flavor. The gap parameters A\, Ai and 
A 3 describe d-s, u-s and u-d quark Cooper pairs, respectively. The following 
pairing schemes have emerged. At asymptotic densities (m s — > or fi — > oo) 
the ground state of QCD with a vanishing strange quark mass is the color- 
flavor locked (CFL) phase (color-flavor locked quark pairing), in which all 
three quark flavors participate symmetrically. The gaps associated with this 
phase are 

A 3 ~A 2 = A 1 =A, (16) 
and the quark condensates of the CFL phase are approximately of the form 

(iilC^-Ae^ef^x, (17) 

with color and flavor indices all running from 1 to 3. Since e a ^ X £f a f b x = 
5J 5j b — <5/ b <5/ one sees that the condensate (17) involves Kronecker delta 
functions that link color and flavor indices. Hence the notion color- flavor lock- 
ing. The CFL phase has been shown to be electrically neutral without any 
need for electrons for a significant range of chemical potentials and strange 
quark masses [112]. If the strange quark mass is heavy enough to be ignored, 
then up and down quarks may pair in the two-flavor superconducting (2SC) 
phase. Other possible condensation patterns are CFL- if [113], CFL-K + and 
CFL-7T '- [114], gCFL (gapless CFL phase) [111], ISC (single-flavor-pairing) 
[111, 115, 116], CSL (color-spin locked phase) [117], and the LOFF (crys- 
talline pairing) [118, 119, 120] phase, depending on m s , (i, and electric charge 
density. Calculations performed for masslcss up and down quarks and a very 
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heavy strange quark mass (m s — > oo) agree that the quarks prefer to pair in 
the two-flavor superconducting (2SC) phase where 

A 3 > , and A 2 = A 1 =0. (18) 

In this case the pairing ansatz (15) reduces to 

(V£C 7 5^ b ><xZU ab e^ 3 . (19) 

Here the resulting condensate picks a color direction (3 or blue in the exam- 
ple (19) above), and creates a gap A at the Fermi surfaces of quarks with the 
other two out of three colors (red and green). The gapless CFL phase (gCFL) 
may prevail over the CFL and 2SC phases at intermediate values of m 2 s /fi 
with gaps given obeying the relation A% > Ai > A\ > 0. For chemical po- 
tentials that are of astrophysical interest, fi < 1000 MeV, the gap is between 
50 and 100 MeV. The order of magnitude of this result agrees with calcu- 
lations based on phenomenological effective interactions [25, 121] as well as 
with perturbativc calculations for \i > 10 GeV [122]. We also note that super- 
conductivity modifies the equation of state at the order of (A/fi) 2 [123, 124], 
which is even for such large gaps only a few percent of the bulk energy. Such 
small effects may be safely neglected in present determinations of models for 
the equation of state of quark-hybrid stars. There has been much recent work 
on how color superconductivity in neutron stars could affect their properties 
[22, 23, 118, 125, 126, 127]. These studies reveal that possible signatures in- 
clude the cooling by neutrino emission, the pattern of the arrival times of 
supernova neutrinos, the evolution of neutron star magnetic fields, rotational 
stellar instabilities, and glitches in rotation frequencies. 

4 Strange Quark Matter 

It is most intriguing that for strange quark matter made of more than a few 
hundred up, down, and strange quarks, the energy of strange quark matter 
may be well below the energy of nuclear matter [15, 16, 17], E/A = 930 MeV, 
which gives rise to new and novel classes of strange matter objects, ranging 
from strangclets at the low baryon-number end to strange stars at the high 
baryon number end [7, 11, 18, 20, 128]. A simple estimate indicates that 
for strange quark matter E/A = ABtt 2 //! 3 , so that bag constants of B = 
57 McV/fm 3 (i.e. B 1 ^ = 145 MeV) and B = 85 MeV/fm 3 (B 1 / 4 = 160 MeV) 
would place the energy per baryon of such matter at E/A = 829 MeV and 
915 MeV, respectively, which correspond obviously to strange quark matter 
which is absolutely bound with respect to nuclear matter [20, 129]. 

4.1 Nuclear crust on strange stars 

Strange quark matter is expected to be a color superconductor which, at 
extremely high densities, should be in the CFL phase [22, 23]. This phase is 
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Fig. 8. Illustration of the eos of strange stars with nuclear crusts (from [130]). 

rigorously electrically neutral with no electrons required [112]. For sufficiently 
large strange quark masses, however, the low density regime of strange quark 
matter is rather expected to form other condensation patterns (e.g. 2SC, CFL- 
K°, CFL-K + , CFL-7r 0, ~) in which electrons are present [22, 23]. The presence 
of electrons causes the formation of an electric dipole layer on the surface of 
strange matter, with huge electric fields on the order of 10 19 V/cm, which 
enables strange quark matter stars to be enveloped in nuclear crusts made of 
ordinary atomic matter [18, 19, 63, 131]. 4 The maximal possible density at the 
base of the crust (inner crust density) is determined by neutron drip, which 
occurs at about 4 x 10 11 g/cm 3 or somewhat below [63]. The eos of such 
a system is shown in Fig. 8. Sequences of compact strange stars with and 
without (bare) nuclear crusts are shown in Fig. 3. Since the nuclear crust is 
gravitationally bound to the quark matter core, the mass-radius relationship 
of strange stars with crusts resembles the one of neutron stars and even that of 
white dwarfs [21]. Bare strange stars obey M oc R 3 because the mass density 
of quark matter is almost constant inside strange stars. 

4.2 Strange dwarfs 

For many years only rather vague tests of the theoretical mass-radius rela- 
tionship of white dwarfs were possible. Recently the quality and quantity of 
observational data on the mass-radius relation of white dwarfs has been re- 
analyzed and profoundly improved by the availability of Hipparcos parallax 

4 Depending on the surface tension of blobs of strange matter and screening effects, 
a heterogeneous crust comprised of blobs of strange quark matter embedded in 
an uniform electron background may exist in the surface region of strange stars 
[132]. This heterogeneous strange star surface would have a negligible electric 
field which would make the existence of an ordinary nuclear crust, which requires 
a very strong electric field, impossible. 
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Fig. 9. Comparison of the theoretical mass-radius relationships of strange dwarfs 
(solid curves) and normal white dwarfs [130]. Radius and mass are in units of Rq 
and Mq, respectively. 



measurements of several white dwarfs [133]. In that work Hipparcos paral- 
laxes were used to deduce luminosity radii for 10 white dwarfs in visual bi- 
naries of common proper-motion systems as well as 11 field white dwarfs. 
Complementary HST observations have been made to better determine the 
spectroscopy for Procyon B [134] and pulsation of G226-29 [135]. Procyon B 
at first appeared as a rather compact star which, however, was later confirmed 
to lie on the normal mass-radius relation of white dwarfs. Stars like Sirius B 
and 40 Erin B, fall nicely on the expected mass-radius relation too. Several 
other stars of this sample (e.g. GD 140, G156-64, EG 21, EG 50, G181-B5B, 
GD 279, WD2007-303, G238-44) however appear to be unusually compact 
and thus could be strange dwarf candidates [136]. The situation is graphically 
summarized in Fig. 9. 

4.3 Surface properties of strange matter 

The electrons surrounding strange quark matter are held to quark matter elec- 
trostatically. Since neither component, electrons and quark matter, is held in 
place gravitationally, the Eddington limit to the luminosity that a static sur- 
face may emit does not apply, and thus the object may have photon luminosi- 
ties much greater than 10 38 erg/s. It was shown by Usov [137] that this value 
may be exceeded by many orders of magnitude by the luminosity of e + e~ 
pairs produced by the Coulomb barrier at the surface of a hot strange star. 
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Fig. 10. Composition of hot (T = 
40 MeV) proto-neutron star matter for 
Y L = 0.3 [142]. 
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Fig. 11. Same as Fig. 10, but for stan- 
dard neutron star matter [142]. 



For a surface temperature of ~ 10 11 K, the luminosity in the outflowing pair 
plasma was calculated to be as high as ~ 3 x 10 51 erg/s. Such an effect may 
be a good observational signature of bare strange stars [137, 138, 139, 140]. 
If the strange star is enveloped by a nuclear crust however, which is gravita- 
tionally bound to the strange star, the surface made up of ordinary atomic 
matter would be subject to the Eddington limit. Hence the photon emissivity 
of such a strange star would be the same as for an ordinary neutron star. If 
quark matter at the stellar surface is in the CFL phase the process of e + e~ 
pair creation at the stellar quark matter surface may be turned off, since cold 
CFL quark matter is electrically neutral so that no electrons are required and 
none are admitted inside CFL quark matter [112]. This may be different for 
the early stages of a hot CFL quark star [141]. 



5 Proto-Neutron Star Matter 

Here we take a brief look at the composition of proto-neutron star matter. 
The composition is determined by the requirements of charge neutrality and 
equilibrium under the weak processes, B% — ► B2 + 1 + H and B2 + I — > Bi + vi, 
where B\ and B2 are baryons, and I is a lepton, either an electron or a muon. 
For standard neutron star matter, where the neutrinos have left the system, 
these two requirements imply that Q — ^\ 1i'nB i + X)z= e u 1i n i = (electric 
charge neutrality) and fiBt — bifi n — qifii (chemical equilibrium), where q^n 
denotes the electric charge density of a given particle, and n^, (ni) is the 
baryon (lepton) number density. The subscript i runs over all the baryons 
considered. The symbol \XBi refers to the chemical potential of baryon i, bi 
is the particle's baryon number, and qi is its charge. The chemical potential 
of the neutron is denoted by /i„. When the neutrinos are trapped, as it is 
the case for proto-neutron star matter, the chemical equilibrium condition is 
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Fig. 12. Representation of an m = 4 perturbation of a rotating neutron star, fi 
denotes the star's rotational frequency, Lu m is the frequency of the counter-rotating 
perturbation [7]. 

altered to [i Bi = &i/x„ - q t (pi - /z„,) and ^ e - fj, Ue = ^ p - /x^, where \x Vl is 
the chemical potential of the neutrino vi. In proto- neutron star matter, the 
electron lepton number Yl = (n e + n Vc )/nB is initially fixed at a value of 
around Yl c = Y e + Y Ve ~ 0.3 — 0.4 as suggested by gravitational collapse cal- 
culations of massive stars. Also, because no muons are present when neutrinos 
are trapped, the constraint Yj, =Y^+ Y VmU = can be imposed. Figures 10 
and 11 show sample compositions of proto-neutron star matter and standard 
neutron star matter (no neutrinos) computed for the relativistic mean-field ap- 
proximation. The presence of the A particle in (proto) neutron star matter at 
finite temperature is striking. This particle is generally absent in cold neutron 
star matter treated in the relativistic mean-field approximation [6, 7, 143]. 

6 Rotational Instabilities 

An absolute limit on rapid rotation is set by the onset of mass shedding from 
the equator of a rotating star. However, rotational instabilities in rotating 
stars, known as gravitational radiation driven instabilities, set a more stringent 
limit on rapid stellar rotation than mass shedding. These instabilities originate 
from counter-rotating surface vibrational modes which at sufficiently high 
rotational star frequencies are dragged forward, as schematically illustrated in 
Fig. 12. In this case gravitational radiation, which inevitably accompanies the 
aspherical transport of matter, does not damp the instability modes but rather 
drives them. Viscosity plays the important role of damping these instabilities 
at a sufficiently reduced rotational frequency such that the viscous damping 
rate and power in gravity waves are comparable. The most critical instability 
modes that are driven unstable by gravitational radiation are f-modes [7, 144] 
and r-modes [145, 146]. Figure 13 shows the stable neutron star frequencies 
if only f-modes were operative in neutron star. One sees that hot as well as 
cold neutron stars can rotate at frequencies close to mass shedding, because 
of the large contributions of shear and bulk viscosity, respectively, for this 
temperature regime. The more recently discovered r-mode instability [145, 
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Fig. 13. Gravitational radiation driven 
f-mode instability suppressed by shear 
and bulk viscosity. (Fig. from [7].) 




Fig. 14. Comparison of f-mode instabil- 
ity with r-mode instability. (Data from 
Refs. [144, 147].) 



148] may change the picture completely, as can be seen from Fig. 14. These 
modes are driven unstable by gravitational radiation over a considerably wider 
range of angular velocities than the f-modes (cf. dashed curve labeled (m = 2) 
r-mode instability). In stars with cores cooler than <~ 10 9 K, on the other 
hand, the r-mode instability may be completely suppressed by the viscosity 
originating from the presence of hyperons in neutron star matter, so that 
stable rotation would be limited by the f-mode instability again [147]. 

Figures 15 and 16 are the counterparts to Figs. 13 and 14 but calculated 
for strange stars made of CFL and 2SC quark matter, respectively [149, 150]. 
The r-mode instability seems to rule out that pulsars are CFL strange stars, if 
the characteristic time scale for viscous damping of r-modes are exponentially 
increased by factors of ~ A/T as calculated in [149]. An energy gap as small 
as A = 1 MeV was assumed. For much larger gaps of A <~ 100 MeV, as 
expected for color superconducting quark matter (see section 3.5), the entire 
diagram would be r-mode unstable. The full curve in Fig. 15 is calculated for 
a strange quark mass of m s = 200 MeV, the dotted curve for m s = 100 MeV. 
The box marks the positions of most low mass X-ray binaries (LMXBs) [151], 
and the crosses denote the most rapidly rotating millisecond pulsars known. 
All strange stars above the curves would spin down on a time scale of hours 
due to the r-mode instability, in complete contradiction to the observation 
of millisecond pulsars and LMXBs, which would rule out CFL quark matter 
in strange stars (see, however, [152]). Figure 16 shows the critical rotation 
frequencies of quark stars as a function of internal stellar temperature for 
2SC quark stars. For such quark stars the situation is less conclusive. Rapid 
spin-down, driven by the r-mode gravitational radiation instability, would 
happen for stars above the curves. 
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Fig. 15. Critical rotation frequen- Fig. 16. Same as Fig. 15, but for 2SC 
cies versus stellar temperature for CFL quark stars [150]. 
strange stars [150]. 



7 Net Electric Fields and Compact Star Structure 

Here we consider the possibility that the electric charge density inside com- 
pact stars (neutron stars, strange stars) is not identically zero. This may be 
the case, for example, for compact stars accreting ionized hydrogen. Another 
example are strange quark stars. They could have electric charge distribu- 
tions on their surfaces that generate electric fields on the order of 10 18 V/cm 
[7, 11, 18, 131] for ordinary quark matter, and 10 19 V/cm [153] if quark matter 
is a color-superconductor. Although the electric field on strange stars exists 
only in a very narrow region of space, it is interesting to study the effects of 
such ultra-high electric fields on the structure of the star. 

It has already been shown that the energy densities of ultra-high elec- 
tric fields can substantially alter the structure (mass-radius relationship) of 
compact stars [154], depending on the strength of the electric field. In con- 
trast to electrically uncharged stars, the energy-momentum tensor of charged 
stars has two key contributions, the usual matter-energy term plus the energy 
density term that originates from the electric field. The latter plays a dual 
role for compact star physics. Firstly, it acts as an additional source of grav- 
ity and, secondly, it introduces Coulomb interactions inside the star. Both 
features can alter the properties of compact stars significantly, as we shall 
demonstrate below. 

We will restrict ourselves to spherically symmetric compact stars. The 
metric of such objects is given by 



ds 2 = e v ^c 2 dt 2 - e x ^dr 2 - r 2 {d9 2 + sin 2 9d<j> 2 ) . 



(20) 
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The energy-momentum tensor consists of the usual perfect fluid term supple- 
mented with the electromagnetic energy-momentum tensor, 



= (p + pc')u K u» + pSS + 



1 

4tt 



F» l F Kl + ±5 K »F kl F k 

47T 



(21) 



where w M is the fluid's four-velocity, p and pc 2 = e are the pressure and energy 
density, respectively, and F^ K satisfies the covariant Maxwell equation, 



(22) 



The quantity J K denotes the four-current which represents the electromagnetic 
sources in the star. For a static spherically symmetric system, the only non- 
zero component of the four-current is J 1 , which implies that the only non- 
vanishing component of F Kft is F 01 . We therefore obtain from Eq. (22) 
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(23) 



which is nothing other than the electric field. This relation can be identified 
as the relativistic version of Gauss' law. In addition we see that the electric 
charge of the system is given by 



Q(r) 



f 4irfe (u 
Jo 



+^' 2 dr. 



(24) 



With the aid of Eq. (24) the energy-momentum tensor of the system can be 
written as 



J- K. 



V 
















Q 2 (r) 



\ 



(25) 



Using the energy-momentum tensor (25), Einstein's field equation leads to 
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(27) 



At this point we define the radial component of the metric g 11 , in analogy to 
the exterior solution of Reissner-Nordstrom, as [155] 



-A (r) = 1 Gm(r) | GQ 2 (r) 
rc 2 r 2 c 4 



(28) 
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From equations (26), (27) and (28), we derive an expression for m(r), which is 
interpreted as the total mass of the star at a radial distance r. This expression 
reads 

dm(r) _ 47rr 2 Q(r) dQ(r) 

— 1 — — — T~ eH 2 7i — ' 

dr c z err dr 

which reveals that, in addition to the standard term originating from the eos 
of the stellar fluid, the electric field energy contributes to the star's total mass 
too. Next, we impose the vanishing of the divergence of the energy-momentum 
tensor, T^ re;/J = 0, which leads to the Tolman-Oppenheimer-Volkoff (TOV) 
equation of electrically charged stars, 
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(30) 



Summarizing the relevant stellar structure equations, we end up with the 
following set of equations: 
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(31) 
(32) 

(33) 
(34) 

(35) 



Equations (31) and (32) arise from Einstein's field equation, Eq. (33) is the 
mass continuity equation, Eq. (34) comes from the Maxwell equations, and 
Eq. (35) is the TOV equation. This system of coupled differential equations 
is subject to the following boundary conditions 



p(0)=p c , e A = 0,Q(0) = 0, m(0) = 0. 



(36) 



In addition to these conditions, one needs to specify the star's central density 
(or, equivalently, the central pressure) for a given equation of state and a given 
electric charge distribution. This will be discussed in more details in the next 
sections. 

As already mentioned at the beginning of this section, strange stars may 
be expected to carry huge electric fields on their surfaces [7, 11, 18, 131, 153]. 
We want to study the effects of such fields on the overall structure of strange 
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Fig. 17. Displacement of electric charges on the surface of a strange star. The 
mathematical form is obtained by superimposing two Gaussian functions. 

stars. To this aim, we model the charge distribution by superimposing two 
Gaussian functions. The first Gaussian is chosen to be positive, representing 
the accumulation of a net positive charge. The second Gaussian, slightly dis- 
placed from the first one, is chosen negative to represent the accumulation of 
a net negative charge. Mathematically, we thus have 

where a is a constant that controls the magnitude of the Gaussians and b the 
widths of the Gaussians. The graphical illustration of Eq. (37) is shown in 
Fig. 17. To obtain a noticeable impact of the electric field on the structure of 
strange stars, one needs to have Gaussians with a width of at least around 
0.05 km. For such widths we find the mass-radius relationships shown in Fig. 
18. The deviations from the mass-radius relationships of uncharged strange 
stars are found to increase with mass, and are largest for the maximum-mass 
star of each stellar sequence. 

The radial distribution of the electric charge over the surface of a strange 
star is particularly interesting. The reason being the occurrence of the metric 
functions in Eq. (24), which defines the star's total net charge. Since the metric 
functions are not symmetric in the radial distance, the charge distribution is 
rendered asymmetric and stars that are strictly electrically charge neutral in 
flat space-time become electrically charged and thus possess non-zero electric 
fields. Figures 19 shows the electric field at the surface of strange stars. Figure 
20 shows the net electric charge at the surface of strange stars. Both plots 
account for the general relativistic charge separation effect. 



2.5 r- 

2 - 

1.5 - 

1 - 

0.5 - 

- 

-0.5 - 

-1 - 

-1.5 - 

-2 - 

-2.5 - 




Neutron Star Interiors and the eos of Superdense Matter 



23 




11.5 11.55 11.6 11.65 11.7 11.75 11.8 

R(km) 



Fig. 18. Mass-radius relationships of electrically charged strange stars. 
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Fig. 19. Electric fields at the surface of strange stars. 



8 Conclusions and Outlook 

It is often stressed that there has never been a more exciting time in the 
overlapping areas of nuclear physics, particle physics and relativistic astro- 
physics than today. This comes at a time where new orbiting observatories 
such as the Hubble Space Telescope (HST), Rossi X-ray Timing Explorer, 
Chandra X-ray satellite, and the X-ray Multi Mirror Mission (XMM New- 
ton) have extended our vision tremendously, allowing us to observe compact 
star phenomena with an unprecedented clarity and angular resolution that 
previously were only imagined. On the Earth, radio telescopes like Arecibo, 
Green Bank, Parkes, VLA, and instruments using adaptive optics and other 
revolutionary techniques have exceeded previous expectations of what can 
be accomplished from the ground. Finally, the gravitational wave detectors 
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Fig. 20. Electric charge on strange quark stars. 



LIGO, LISA, VIRGO, and Geo-600 are opening up a window for the de- 
tection of gravitational waves emitted from compact stellar objects such as 
neutron stars and black holes. This unprecedented situation is providing us 
with key information on neutron stars, which contain cold and ultra-dense 
baryonic matter permanently in their cores. As discussed in this paper, a key 
role in neutron star physics is played by strangeness. It alters the masses, radii, 
moment of inertia, frame dragging of local inertial frames, cooling behavior, 
and surface composition of neutron stars. Other important observables influ- 
enced by strangeness may be the spin evolution of isolated neutron stars and 
neutron stars in low-mass x-ray binaries. All told, these observables play a key 
role for the exploration of the phase diagram of dense nuclear matter at high 
baryon number density but low temperature [184], which is not accessible to 
relativistic heavy ion collision experiments. 

Obviously, our understanding of neutron stars has changed dramatically 
since their first discovery some 40 years ago. In what follows, I briefly sum- 
marize what we have learned about the internal structure of these fascinating 
object since their discovery. I will address some of the most important open 
questions regarding the composition of neutron star matter and its associated 
equation of state, and will mention new tools, telescopes, observations, and 
calculations that are needed to answer these questions: 

• There is no clear picture yet as to what kind of matter exists in the cores of 
neutron stars. They may contain significant hyperon populations, boson 
condensates, a mixed phase of quarks and hadrons, and/or pure quark 
matter made of unconfined up, down, and strange quarks. 

• Pure neutron matter constitutes an excited state relative to many-baryon 
matter and, therefore, will quickly transform via weak reactions to such 
matter. 
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Table 1. Past, present, and future search experiments for strange quark matter 
[11]. 
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Search for strangelets in terrestrial matter: 


[170] 


Tracks in ancient mica 


[171, 172] 


Rutherford backscattering 


[173, 174] 


Search for strangelets at accelerators: 




Strangelet searches E858, E864, E878, E882-B, 


[175, 176, 177] 


E896-A, E886 




H-dibaryon search 


[178, 179] 


Pb+Pb collisions 


[180, 181, 182, 183] 



a AMS: Alpha Magnetic Spectrometer (scheduled for 2005-2008). 

b CRASH: Cosmic Ray And Strange Hadronic matter. 

c ECCO: Extremely-heavy Cosmic-ray Composition Observer. 

d 1MB: Irvine Michigan Brookhaven proton-decay detector (1980-1991). 

° JACEE: Japanese- American Cooperative Emulsion Chamber Experiment. 

f MACRO: Monopole, Astrophysics and Cosmic Ray Observatory (1989-2000). 



• Neutron stars made up of pure, interacting neutron matter cannot rotate 
as rapidly as the very recently discovered pulsars PSR J1748-2446ad, which 
spins at 716 Hz. The equation of state of such matter, therforc, imposes 
an upper bound on the equation of state of neutron star matter that is 
tighter than the usual P — e constraint (see Fig. 2). 

• Charm quarks do not play a role for neutron star physics, since they be- 
come populated at densities which are around 100 times greater than the 
densities encountered in the cores of neutron stars. While hydrostatically 
stable, "charm" stars are unstable against radial oscillations and, thus, 
cannot exist stably in the universe [131]. 

• Multi-quark states like the H-particle appear to make neutron stars un- 
stable. 

• Significant populations of A's are predicted by relativistic Brueckncr- 
Hartree-Fock calculations, but not by standard mean-field calculations 
which do not account for dynamical correlations among baryons computed 
from the relativistic T-matrix equation. 

• The finite temperatures of proto neutron stars favors the population of 
Zi's already at the mean- field level. 
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• The r- modes are of key interest for several reasons: 1. they may explain 
why young neutron stars spin slowly 2. why rapidly accreting neutron 
stars (LMXB) spin slowly and within a narrow band, and 3. they may 
produce gravitational waves detectable by LIGO. Knowing the bulk vis- 
cosity originating from processes like n + n— > p + + E~ and the superfluid 
critical temperature of S~ , both are poorly understood at present, will be 
key. 

• The loss of pressure resulting from the appearance of additional hadronic 
degrees of freedom at high densities reduces the (maximum) mass of neu- 
tron stars. This feature may serve as a key criteria to distinguish between, 
and eliminate certain, classes of equations of state [6, 7, 185]. 

• Heavy neutron stars, with masses of around two solar masses, do not au- 
tomatically rule out the presence of hypcrons or quarks in the cores of 
neutron stars [186]. 

• Depending on the densities reached in the cores of neutron stars, both 
Schroedinger-based models as well as relativistic ficld-thcorctical models 
may be applicable to neutron star studies. 

• The density dependence of the coupling constants of particles in ultra- 
dense neutron star matter needs be taken into account in stellar structure 
calculations. Density dependent relativistic field theories are being devel- 
oped which account for this feature.. 

• The models used to study the quark-hadron phase transition in the cores 
of neutron stars are extremely phenomenological and require considerable 
improvements. 

• If quark matter exists in the cores of neutron stars, it will be a color 
superconductor whose complex condensation pattern is likely to change 
with density inside the star. The exploration of the numerous astrophysical 
facets of (color superconducting) quark matter is therefore of uppermost 
importance. What are the signatures of color superconducting quark mat- 
ter in neutron stars? So far it has mostly been demonstrated that color 
superconductivity is compatible with observed neutron star properties. 

• A two-step quark-hadron phase transition (1. from nuclear matter to reg- 
ular quark matter, 2. from regular quark matter to color superconducting 
quark matter) may explain long quiescent gamma-ray bursts due to the 
two phase transitions involved. 

• Are there isolated pulsar that are spinning up? Such a (backbending) phe- 
nomenon could be caused by a strong first-order-like quark-hadron phase 
transitions in the core of a neutron star [105, 187, 188]. 

• Was the mass of the neutron star created in SN 1987A around 1.5M Q ? 
And did SN 1987A go into a black hole or not? If the answer to both ques- 
tions were yes, a serious conflict with the observation of heavy neutron 
stars would arise. On the other hand, it could also indicate the existence 
generically different classes of "neutron" stars with very different maxi- 
mum masses. 
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Sources known to increase the masses of neutron stars are differential ro- 
tation, magnetic fields, and electric fields. Some of these sources are more 
effective (and plausible) than others though. 

Nuclear processes in non-equilibrium nuclear crusts (e.g. pycnonuclear re- 
actions) and/or cores (heating caused by changes in the composition) of 
neutron stars can alter the thermal evolution of such stars significantly. 
We are just beginning to study these processes in greater detail. 
What is the shell structure for very neutron rich nuclei in the crusts of 
neutron stars? 

Do N=50 and N=82 remain magic numbers? Such questions will be ad- 
dressed at GS1 (Darmstadt) and R1KEN. 

Are there pulsars that rotate below one millisecond? Such objects may be 
composed of absolutely stable strange quark matter instead of purely grav- 
itationally bound hadronic matter. Experimental physicists have searched 
unsuccessfully for stable or quasistable strange matter systems over the 
past two decades. These searches fall in three main categories: (a) searches 
for strange matter (strange nuggets or strangelets) in cosmic rays, (b) 
searches for strange matter in samples of ordinary matter, and (c) at- 
tempts to produce strange matter at accelerators. An overview of these 
search experiments is given in table 1. 

Strange stars may be enveloped in a crust. There is a critical surface 
tension below which the quark star surfaces will fragment into a crystalline 
crust made of charged strangelets immersed in an electron gas [128, 132] 
If bare, the quark star surface will have peculiar properties which distin- 
guishes a quark star from a neutron star [137, 138, 189, 190]. 
A very high-luminosity flare took place in the Large Magellanic Cloud 
(LMC), some 55 kpc away, on 5 March 1979. Another giant flare was 
observed on 27 August 1998 from SGR 1900+14. The inferred peak lumi- 
nosities for both events is ~ 10 7 times the Eddington limit for a solar mass 
object, and the rise time is very much smaller than the time needed to drop 
~ 10 25 g (about 10 _8 M Q ) of normal material onto a neutron star. Alcock et 
al. [18] suggested a detailed model for the 5 March 1979 event burst which 
involves the particular properties of strange matter (see also [190, 191]). 
The model assumes that a lump of strange matter of <~ 10 -8 M Q fell onto 
a rotating strange star. Since the lump is entirely made up of self-bound 
high-density matter, there would be only little tidal distortion of the lump, 
and so the duration of the impact can be very short, around ~ 10~ 6 s, 
which would explain the observed rapid onset of the gamma ray flash. 
The light curves expected for such giant bursts [137, 138, 139, 140] should 
posses characteristic features that are well within the capabilities of ESA's 
INTErnational Gamma- Ray Astrophysics Laboratory (INTEGRAL [192]) 
launched by the European Space Agency in October of 2002. 
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